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Using the proton ordering model modified by taking into acconnt the dependence 
of the dipole moments on the order parameter [Vdovych et al, 2014], we explore the 
piezocaloric and mnlticaloric effects in the KH2PO4 type ferroelectrics, cansed by 
the shear stress and longitndinal electric field E^. The mnlticaloric effect is 
shown to be stronger than either electrocaloric or piezocaloric effects, especially at 
temperatnres far from the phase transition. 
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1 Introduction 

The piezocaloric effect is a change of the sample temperature AT at adiabatic chang¬ 
ing of the applied to it mechanical stress. This effect, just like the electrocaloric 
effect, is explored because of its potential application in environment-friendly and 
compact solid-state cooling devices. 

The largest piezocaloric effect so far has been observed in shape memory alloys, 
undergoing a martensitic transformation. Thus, in NiTi it was obtained that AT = 
17 K under the tensile stress of cr ~ 650 MPa [1] and AT = 40 K at a ~ 800 MPa 
[2]. In Cu68.i3Zni5,74Ali6,i3 the value of AT = 15 K was reached [3]. 

Ferroelastics are promising materials too. The combined electrocaloric and 
piezocaloric effects (a mnlticaloric effect) have been studied [1] by the ab initio 
method in PbTiOs crystals, which are both ferroelectric and ferroelastic. The cal¬ 
culated temperature change in large fields and stresses exceeds 30 K. 

The hydrogen bonded ferroelectrics of the KH2PO4 (KDP) family are also fer¬ 
roelastic. In the ferroelectric phase the spontaneous polarization P3 is accompanied 
by the spontaneous shear strain e^. So far, neither the piezocaloric nor mnlticaloric 
effcts in these crystal have been explored yet. 

Influence of the shear stress Ug and electric field T3 on polarization, dielectric 
permittivity, piezomodules, elastic constants of the KH2PO4 family crystals was de¬ 
scribed in [5l [6], [71 [8] using the proton ordering model modified by taking into account 
a piezoelectric coupling to the shear strain Sq. Shifts of the transition temperature 
and eventual smearing out of the first order transition by these fields conjugate to 
the order parameter were also described. However, these theories suffered from the 
inner logical contradiction, associated with the required invoking of two different 
values of the effective dipole moments for the paraelectric and ferroelectric phases 
Pi- Hence, while no physical characteristic of the crystal should exhibit any 
discontinuity in the fields above the critical one, there was no smooth transition 
between the values of model parameters, set to be different for the two phases. 
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This contradiction has been removed in (Qj [10] by taking into account the assumed 
dependence of the effective dipole moment on the order parameter. The term 

has been introduced to the Hamiltonian, equivalent to a term proportional to 
in a phenomenological thermodynamic potential. Such a modihcation allowed us 
POTO] to quantitatively and consistently describe the behavior of the physical char¬ 
acteristics of the KH 2 PO 4 family crystals in presence of the electric held includ¬ 
ing the electrocaloric effect. 

In this paper we use the model [10] to explore the piezocaloric and multicaloric 
effects in the KH 2 PO 4 type crystals as a reaction to the shear stress ag. 


2 Calculations 

No modihcation to the model developed in [10] is required to study the inhuence of 
(Tg stress. All the derived expressions for the entropy, polarization, and shear strain 
eg remain valid and can be found in P [TU] . 

From the entropy we obtain the molar specihc of the proton subsystem = 
T{dS/dT)„. The total specihc heat is considered to be the sum of AC and of a 
some regular background contribution, mostly from the lattice of heavy ions 

C = AC + Creg ( 2 ) 

The regular contribution near T^. is approximated by the linear dependence 

Ce 3 = C 0 + C'i(T-Tc), 

yielding a good agreement with experimental data PEQ]. The corresponding regular 
contribution to the entropy near is then 

Sreg = j ^dT = (Cg - CCc) ln(T) + cc + const 

Hence, the total entropy as a function of the held and stress cxg is 

Stotal{TC.(T&) = S + Sreg. ( 3 ) 

Solving (P) with respect to temperature at StotaiiT, E, cxg) = const and two values of 
the held and stress, one can calculate the electrocaloric, piezocaloric, or combined 
multicaloric temperature shifts (as seen in hg. [ 1 ]) 

AT = nStotai, E{2), ag(2)) - T{Stotai, E{1), ag(l)). ( 4 ) 

Numerical calculations are performed for the K(Hi_ 3 ;Da;) 2 P 04 crystals (at a; = 0 
undergoing the hrst order phase transition close to the second order transition) and 
K(Hi_a;Da,) 2 As 04 (with a pronounced hrst order transition at x = 1). For this we 
use the same values of the model parameters that were obtained in P [ 10 ] . 
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Behavior of molar entropy (fig. [T]), polarization and shear strain Eq (hg. [ 2 ]) is typ¬ 
ical for the systems with the hrst order phase transitions in external helds conjugate 
to the order parameter. The stress dg induces non-zero strain Eq and polarization P 3 
above the transition temperature. At low stresses both have jumps at Tc, indicating 
the hrst order phase transitions. The jumps decrease with increasing stress, and Tc 
increases. Above the certain critical stress CTg (the corresponding critical tempera¬ 
ture is T*) the phase transition is smeared out. The critical stress and temperature 
are = 92 kPa, T* = 122.246 K for a; = 0 and = 6.2 MPa, T* = 212.585 K for 
X = 0.89. 




Figure 1: The temperature dependence of molar entropy of KH2PO4 (left) and 
KD 2 ASO 4 (right) at different values of the stress ag. 



Figure 2: The temperature dependence of polarization and shear strain Eq of 
KH 2 PO 4 at different values of the stress cxg. 

The calculated dependence of the piezocaloric temperature change AT of the 
studied crystals on the adiabatically applied stress erg is shown in hgs. [3] and HI 

At small CTg the piezocaloric temperature change increases with the stress linearly 
at temperatures below T° (curves 1, 1’) and quadratically above T° (curves 3, 3’). 
Here T° denotes the transition temperature in absence of external helds. At T = 
TO and CTg < CTg the temperature change is equal to the shift of the transition 
temperature, which is roughly proportional to the stress (see hg. [1]). It is well seen 
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Figure 3: The stress dependence of the piezocaloric temperature change of 
K(Hi_j.Dj.) 2 P 04 crystals at a: = 0.0 (solid lines) and x = 0.89 (dashed lines) at 
T - TO = -5 K - 1, T = TO - 2, T - TO = 5 K - 3. 



Figure 4: The same of KH 2 ASO 4 (solid lines) and KD 2 ASO 4 (dashed lines) atT — 
TO = -2 K ~ 1, T = TO - 2, T - TO = 2 K - 3. 

at (Jg < 40 MPa in the KD 2 ASO 4 crystal, undergoing the most pronounced hrst 
order phase transition (£g. 01 curve 2’). In KH 2 PO 4 , where the phase transition in 
zero helds is close to the second order, at T = and small stresses ag < 10 MPa the 
temperature change follows the law AT ~ . At large stresses the piezocaloric 

stress above TO is larger than below it, and the deviation from the linear or quadratic 
behavior is observed; at ug ^ 200 MPa the saturation is reached. Unfortunately, 
no experimental data for AT is available. 

As follows from the temperature dependence of AT (£g. |5]), the largest 
piezocaloric temperature change is expected just above Tc, where it can exceed 
5 K. 

Simultaneous application of the electric held T 3 and stress cjg yields the mul¬ 
ticaloric effect (£g. [ 6 ]). It is stronger than either electrocaloric or piezocaloric effect, 
and at temperatures far from T° the multicaloric temperature shift is even larger 
than the simple sum of the respective electro- and piezocaloric shifts. 
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Figure 5: The temperature dependence of the piezocaloric temperature change of 
K(Hi_a,Dj,) 2 P 04 at X = 0.0 (left) and x = 0.89 (right) at different values of the 
stress (Tg. 



Figure 6 : The temperature dependence of the temperature change of the KH 2 PO 4 
crystal under the electric held m (curves 1 ), stress Ug (curves 2 ), the held and 
stress simultaneously (curves 3). 

3 Conclusions 

Modihcation of the proton ordering model, suggested earlier in order to describe the 
electrocaloric ehect in the KH 2 PO 4 family crystals, is used to explore the piezocaloric 
and multicaloric ehects. They are found to be qualitatively and quantitatively sim¬ 
ilar to the electrocaloric ehect. The theory predicts that the temperature change 
can exceed 5 K at high stresses (obviously well above the sample strength) and can 
be increased further by a simultaneous application of both electric held and shear 
stress CTg. Experimental measurements of AT are required. 
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